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A single qubit driven by an appropriate sequence of control pulses can serve as a spectrometer 
of local noise affecting its energy splitting. We show that by driving and observing two spatially 
separated qubits, it is possible to reconstruct the spectrum of cross-correlations of noises acting 
at various locations. When the qubits are driven by the same sequence of pulses, real part of 
cross-correlation spectrum can be reconstructed, while applying two distinct sequence to the two 
qubits allows for reconstruction of imaginary part of this spectrum. The latter quantity contains 
information on either causal correlations between environmental dynamics at distinct locations, or 
on the occurrence of propagation of noisy signals through the environment. We illustrate the former 
case by modeling the noise spectroscopy protocol for qubits coupled to correlated two-level systems. 
While entanglement between the qubits is not necessary, its presence enhances the signal from which 
the spectroscopic information is reconstructed. 


A qubit interacting with its environment experiences 
decoherence that limits the timescale on which it 
can be used for quantum information processing pur¬ 
poses. The time dependence of coherence decay is de¬ 
termined by dynamics of the environmental degrees of 
freedom coupled to the qubit. When the environment is 
well characterized, decoherence is simply a nuisance. On 
the other hand, when the dominant source of decoher¬ 
ence is unknown, measurements of qubit’s coherence de¬ 
cay can be used to obtain substantial information about 
environmental fluctuations. 

Here we focus on the case in which the environment 
is a source of classical noise that affects the energy split¬ 
ting of the qubits, i.e. it leads to pure dephasing. Driving 
the qubit with a sequence of dynamical decoupli ng ( DD) 
pulses [Ml not only slows down decoherence (lll - [l^ . 
but for an appropriately chosen (essentially periodic) se¬ 
quence of n pulses, the magnitude of qubit’s coherence 
at a given time t is proportional to spectral density of 
noise, S'(w), evaluated at w = nn/t [isl - ll^ . Since the 
application of DD pulses to the qubit translates to mod¬ 
ulating the phase noise with a periodic piece-wise con¬ 
stant function of alternating sign, this can be most e asily 
understood as noise filtering by a lock-in mechanism [l8l| . 

The efficacy of DD-based environmental noise spec¬ 
troscopy (DDENS) with a single qubit was shown in 
many experiments on various kinds of qubits, including 
those based on trapped ions [Hill, su perc onducting 
circuits [l^ , semiconductor quantum dots |2Cll . [2lj , phos¬ 
phorous donors in silicon [22|, and NV centers in dia¬ 
mond [1^ [13 . It is crucial to note that in the case of 
solid-state based qubits, each qubit is interacting with a 
specific nanoscale environment, the exact properties of 
which vary from qubit to qubit. Furthermore, qubits 
tightly localized in a nanostructure (e.g. NV centers lo¬ 
cated close to a surface of a diamond nanocrystal), can 
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be brought into contact with various environments, al¬ 
lowing for nanoscale resolution sensing of environmental 
fluctuations 1^1^. 

It is natural to ask if using more than one qubit brings 
new features to DDENS (similarly as in the case of quan¬ 
tum metrology in which using multi-qubit entan¬ 

gled states enhances the signal sensing precision). Mo¬ 
tivated by recent experimental success of single-qubit 
DDENS, we consider its simplest (and thus currently 
realistic experimentally) extension to the case of two 
qubits, each coupled to distinct noise. This way we can 
gain access to cross-correlation of the two noises, which 
enhances our knowledge of the environmental dynamics 
when qubits are close enough to each other so that their 
respective environments are nontrivially correlated. 

We show how from measurements of an appropriate 
two-qubit coherence (e.g. by two-qubit tomography), one 
can reconstruct the spectrum of cross-correlations of the 
noises felt by the two qubits. The real part of this so- 
called cross-power [2^ quantifies the degree to which the 
noise at a given frequency is common for the two qubits. 
This quantitity can be reconstructed using the same DD 
sequence applied to the two qubits. However, with an 
appropriate choice of two distinct sequences, one gains 
access to an imaginary part of the cross-spectrum. The 
latter quantity contains information on time-asymmetric 
correlations between the two noises. These correlations 
can arise either due to causal relation between the two 
noises (which is a signature of interaction among the en¬ 
vironmental degrees of freedom responsible for the two 
noises), or due to the propagation of a noisy signal from 
one part of the environment to another (e.g. due to pres¬ 
ence of common cause affecting one environment with a 
delay with respect to another). While the relevant mea¬ 
surement signal is maximal when the two qubits are ini¬ 
tially entangled, for an appropriate separable state the 
amplitude of the coherence which needs to be recon¬ 
structed can be only two times smaller. In the latter case 
the spectroscopy scheme presented here should not be 
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qualitatively harder to implement than the single-qubit 
procedure - we only need separate but synchronized read¬ 
out of the qubits (for state tomography) and the possi¬ 
bility of applying distinct sequences of tt pulses to each 
of the qubits. 

We consider pure dephasing of qubits 1 and 2: 

^ (f 2 «+Ca(t))<Ti“V 2 , ( 1 ) 

a=l,2 

where fia is the splitting of qubit a, and ^a(t) is the clas¬ 
sical noise affecting this splitting (with average {^a{t))=0 
without any loss of generality). The matrix of two- 
point noise correlation functions is given by Cap(t) = 
(^a(t)^^(O)), and the corresponding matrix of spectral 
densities is Safi{oj) = e*‘^‘C'a/ 3 (t)dt , where Saaiijj) 

is the self-power (first spectral density) of noise ^a(t). 
Since Capit) =Cpa{—t), the cross-correlation spectrum, 
Sapii^) with l3, can be written as 

Sc.p (w) = ^^3(w) + ^Sipico) = S*0^ico) , (2) 

where S^p(uj) and S^p{uj) are real functions (even and 
odd in w, respectively), which characterize the even (odd) 
in t parts of Capit). 

Let us illustrate the physical meaning of Sl^p{uj) and 
S^p{uj) with simple examples. When noises ^a{t) are 
caused by multiple sources (e.g. fluctuating spins or 
charges) located in vicinity of the qubits, we have ^ait) = 
where is the coupling of qubit a to 
the fc-th source of noise. If the sources are uncorre¬ 
lated, {Vk{t)r]i{0)) = SkiCkiit), then Ci 2 {t) = C 2 i{t) = 
'^^ 2 '’Ckk and Sl^p{uj) = 0 , while S^p{uj) is a sum 
of spectra of r]k noises weighed by the product of cou¬ 
pling constants to the two qubits, i.e. it contains infor¬ 
mation on the amount of common noise felt by the qubits. 
Nonzero S^p{uj) will appear when r]k noises become 
causally correlated, i.e. when there is interaction between 
the entities the dynamics of which is the source of noise 
[13 )• For example, when noise ry at given time influences 
the subsequent fluctuations of r]k (but not vice versa), we 
have Cki (t) ^ 0 while cik {t) = 0, and C 12 (t) ^ C 21 (t) if both 
qubits are coupled to r]k and iji . Another simple situation 
in which a finite S^p{uj) appears is that of signal prop¬ 
agation: when the qubits feels the same noise, but with 
distinct delay times i.e. ^a{t) = Va^o{t — we 
obtain 5 'i 2 (a;) = niW 2 e*‘^^* 5 'o(w), where At = t^^^ — t^‘^\ 
and 5*0 (w) is the self-spectrum of fo noise. 

We assume that the two qubits are initialized in state 
/3(0). Then, during the evolution time T each of them 
is subjected to a sequence of ideal tt pulses about its x 
or y axes. We assume that each of the qubits can be 
separately addressed, so that it is possible to apply two 
distinct sequences of pulses, each parametrized by a time- 
domain filter function [ 3 , |3l| 


which is alters between ±1 for 0 ^ t < T and is zero oth¬ 
erwise. The sign changes at the times (fc = l,...,n„) 

at which subsequent pulses are applied to a given qubit 
(note that we define 5o = 0 and 6n+i = 1). 

Even though we do not require ^a(t) to be Gaussian, 
we will treat them as such from now. The justification 
is that in most cases the non-Gaussian noise becomes 
effectively Gaussian upon filtering Biillll (one ex¬ 
ception being the case of quadratic coupling to a very 
low-frequency noise [13 )• The averaging over the noise 
realizations of density matrix p{T) evolving due to Hamil¬ 
tonian from Eq. © and applied pulses, is then a straight¬ 
forward generalization of the calculation of single-qubit 
dephasing due to Gaussian noise 0, mi . The matrix el¬ 
ements of p(T), written in a standard basis of products 
of eigenstates of 0 - 1 “^ are given by 


Pcri(y2,cr'i<y'2 (T) 

= ((-l)"^ai, (-1)-V2|p(0)|(-1)-Vi, (-1)-V')x 

rT 


X ( exp 


“ 2 ^ -'^a) / fa{t')^a{t')dt' 

a = l ,2 -^0 


(4) 


where a a = ±1- We notice that in order for the result 
to be influenced by both and ^2 we need to focus on 
coherences with cr(, = — ctq. Galculating the standard 
Gaussian average for those cases we get 

-ai-a 2 (T) oc exp(-xii - X 22 - 2cri 0 - 2 X 12 ) , (5) 

where the quantities are well-known from the case 
of the single qubit 0, : 

Xo.a = ll^ ( 6 ) 


where fa(pj) is the Eourier transform of fait). 

The key feature of /a(w) functions (or rather \fa{uj)\‘^) 
that enables noise spectroscopy is that for periodic se¬ 
quence of many pulses (n Id they act like narrow- 
band frequency filters Bin [ 3 . In what follows, we 
shall consider two, particular sequences, namely: Garr- 
Purcell (CP) sequence f^p, defined by Sk = {k — ^)/n, 
and the periodic dynamical decoupling (PDD) sequence 
fpDD’ defined by 6k = k/{n + I). Eor odd n, those se¬ 
quences are related in a following way 


7n^T ( \ 

/pDDV^) ~ 


T 


n + 1 


-sine 


^n+l.T / \ —i \ 

fcp (‘^) = e ^^/pbD 


r luT \ Sin(^) 

V2(n+l);cos(2^) 

(w) -I- 0{n~^). 


(7) 

( 8 ) 


Hence, the filtering term that appears under the integral 
in Eq. dH]) can be approximated by (see Fig. [T]) 


ut) = ^(-i)'=0(4?iT - mt - 


(3) 


I/pddMI^-£^ 


4r 




m=±(l,3,...) 


—I a; — m 

'm ^ 


(n + l)7r 

T 


(9) 
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The spectroscopy of 5^2 (w) requires applying two dis¬ 
tinct sequences of pulses to the two qubits. In order for 
this quantity to contribute to integral (TO, the filtering 
term /i(—a;)/2(w) has to be odd in ui. A possible con¬ 
figuration which we propose here is to set /i = /pf^ 
and /2 = /cp^’^ (or rifce versa). Upon inspection of 
Eqs. (O and dS]), we see that this filtering function has 
the desired property due to the relative phase difference 
between FDD and CP sequences 


/pb^D(w)/Sp^’’^(-w) = l/pm3(w)P = 

4T A (n-|-l)7r' 




( 13 ) 


FIG. 1. Illustration of the filtering term |/pi^(a))P (thick, 
black). Function sin^(^)/cos^( ) (bine) has a 

form of a series of very narrow peaks of height located at 
frequencies m(n-|-l)7r/r (m = ±1, ±3,...). Filtering term is a 
product of this function and sinc^)^^^^) (dashed, red) which 
dampens the amplitude of m-peak by a factor of ^jijp'rr?). 


Here ^{z) is a very narrow function (with the width of 
order n~^) centered at z = 0 with height of order n, that 
satisfies A(2;)dz = 27 r. Effectively, the filter (j^ be¬ 
haves as a series of Dirac delta-like functions at frequen¬ 
cies uJm = m{n + 1 )tt/T ( m = ± 1 , ± 3 , ± 5 ,... ), with 
each term weighted by m~^. The most crude approxi¬ 
mation is to keep only the first two terms in series 
that correspond to m = ±1 peaks. Setting fa = /pf^ 
or fa = /cp^’^ s-iid employing this approximate form in 
Eq. dH]) gives us 


Xaa 


^ f (rz_+j>\ 
^2 T ) 


( 10 ) 


which is the basis of the curre ntly most popular single¬ 
qubit DDENS method [H, EE [ij El [13 

In the case of two-qubit coherence decay, the new con¬ 
tribution in Eq. (0) is 

1 poo 1 

X12 = 2 y -S'l2(w)/l(-w)/2(w)^ , (11) 


which encapsulates the influence of noise cross¬ 
correlations on dephasing. We focus now on this term 
exclusively. 

When the same sequence of pulses is applied to both 
qubits, i.e. when /i(-w)/2(w) = /(-a;)/(w) = |/(w)p, 
the hltering term in Eq (HD is even in lo and we ob¬ 
tain an expression for Xi2(t) which is exactly analogous 
to Eq. (ED, only with replacing Saa{iz). Conse¬ 

quently, for CP or PDD sequence discussed above the 
spectroscopic formula for the real part of cross-power is 


{/i-/2-/pdd’^Cp} : Xl2 


A sS 


• (12) 


Keeping only the first two terms corresponding to the 
largest peaks we get 



^n,T 
' /pDD 

rn+l,T 

JCP 


Xl2 


_^qi ( iri +l)Tr \ 
7r2 T ) 


, ( 14 ) 


The above formula, together with Eq. dl2D, are the main 
results of this paper. Utilizing those relations it should 
be possible to perform spectroscopy of cross-correlations 
of two noises felt by two spatially separated (possibly 
distant) qubits. We envision the practical implementa¬ 
tion of the spectroscopy of the cross-power in a following 
procedure. 

1 . Firstly, the decay of the two-qubit coherence, 
Pa^a^-ai-aAT) = exp[-r(T, n)] accompanied by 
pulse sequences fi = /pf^ and /2 = /cp should 
be measured (see Eq. (ED). For large number of 
pulses n this would result in 


r('r, n)=xii+x22+'^(Ti(T2Si2{ ^^ 


It is important to note that while the initial am¬ 
plitude of the coherence is maximal for an entan¬ 
gled two-qubit state (i.e. when a singlet state, or 
|'I'_) Bell state, is created, we have p+i-i .-i+i(0) = 
1 /2), for a separable state such as |-|-cc)^ 0 \+x)2 

(where |-|-x)^ is an eigenstate of (7i“^) we have 
PCTicr2,-o-i-CT2( 0 ) = 1 / 4 - Entanglement is thus help¬ 
ful, since it provides a larger signal, but it is not 
necessary, and one should weigh its benefits against 
effort which typically accompanies creation of two- 
qubit entanglement in a solid-state setting . 

We conjecture that with more than two qubits em¬ 
ployed, the benefits of using an entangled state 
over a separable one will be more pronounced, see 

e.g. M- 


2. In the second step, the measured exponent r(T, n) 
should then be corrected for the presence of single¬ 
qubit dephasing Xaa- This can be accomplished 
in two ways. First one is to characterize the 
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self-power spectra Saai^^) by separately employing 
single-qubit DDENS method and then calculating 
the contribution to T(T,n) according to Eq. ([6|). 
Alternatively, bare Xaa can be extracted from the 
measurement of the decay of an appropriate co¬ 
herence, provided the consistent choice of pulse 
sequence and experiment time. For example, ac¬ 
cording to Eq. (U]), measurements of paia2,-<Ti<T2{T) 
would yield xiij etc.. 

Repeating the procedure with different choices of n 
and T provides access to cross-power at wide range of 
frequencies. The real part of the cross-power can be ex¬ 
tracted in a analogous scheme, with the only difference 
that both qubits have to be treated with the same pulse 
sequence. 

We illustrate now the performace of the above spec¬ 
troscopy scheme by using it to reconstruct the noise 
spectra from a simulation of two nontrivially correlated 
noises. The example that we use is inspired by physics 
of superconducting qubits, which are often strongly af¬ 
fected by random telegraph noise (RTN) generated by 
two-level systems (TLSes) ubiquitous in condensed mat¬ 
ter environment. Experimental signatures of interactions 
between the TLSes have been a subject of recently in¬ 
creased attention [ssl - li^ . Here we focus on a simple 
model in which one TLS (labelled A) strongly affects the 
other TLS (labelled B). The levels of these systems are 
labeled 0 and 1, and the 0—>-1 (1—>0) transition rates of 
a given TLS are y-j- (74,). We will assume that TLS A has 
a switching rate XA = , i-e. its energy splitting 

Ayi is much smaller than ksT, and it is weakly affected 
by the state of TLS B. For the latter one we assume that 
its switching rates 7^ and 7^ do depend on the state of 
A given by a = 0 , 1 . 

For simplicity we assume that TLS A (B) is coupled 
only to qubit 1 ( 2 ). This immediately excludes the sim¬ 
plest reason for nonzero cross-correlation, which is due to 
presence of common noise component. The noises felt 
by the qubits are then simply proportional to the two 
RTN signals. All the Cap{t) correlation functions can 
be calculated exactly (see appendix [B]). While the auto¬ 
correlation of is simply given by the standard RTN 
2 

formula, C'ii(t) = the cross-correlation Ci2(t) 

is nontrivial: 

= ( 16 ) 

where p®*), are the average occupations of the states of 
the two TLS labeled by a,h — 0 , 1 . We see that when 
the average state of B TLS depends on state of A (which 
occurs when 7it/7u?^7ot/7o4.), we have C'i2(t)7^0. The 
remaining correlation functions, C2i(t), and C22(i) are 
given by rather lengthy expressions (see appendix El). 
Since C'i2(t) 7^ C2i(t), we have nonzero 3(2(1^), which 
results from interaction between the two TLS. 

We have generated the noise from the two TLS using 
standard methods, and used it to simulate the whole pro¬ 


cedure of DDENS described above. In addition we also 
checked for corrections due to filter peaks at higher fre¬ 
quencies [1^ (see Eq. (j9] and appendix lAl) and Fig. [T|). 
The results are shown in Fig. [21 where one can see accu¬ 
rate reconstruction of nontrivial shape of 5(2 ((^)- Note 
that while the noise statistics is in fact non-Gaussian, 
under the application of many pulses the noisy phase 
becomes effectively a Gaussian variable [3, [l^ , and 

the theory derived above under assumption of Gaussian 
noise statistics very well describes the coherence decay. 
The negligence of non-Gaussian contributions introduces 
only small systematic error that overstates the value of 
cross-power, especially for low frequencies. 
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FIG. 2. The reconstruction of imaginary and real (inset) 
parts of cross-power in a numerical experiment. The values 
of cross-power have been extracted from data obtained from 
measurement of with (blue, filled circles) and with¬ 

out (red, empty circles) corrections from higher-frequency 
peaks in filtering term (see appendices O and EJ- The re¬ 
sults are compared with the exact function (solid black line). 
The parameters of the coupled telegraph noises were set to: 
7^ = 7u = 7^. = 7n = 7 a = 1 and vi = V 2 ^ 0.1 . The 

number of pulses n = 37 have been kept fixed and the total 
time T was manipulated in order to sweep the wide range of 
frequencies. 

In conclusion, we proposed a realistic method for ex¬ 
tending the experimentally successful single-qubit noise 
spectroscopy techniques to the case of two qubits, each 
affected by classical phase noise. With a proper choice 
of pulse sequence applied to each of the qubits, both real 
and imaginary parts of the cross-correlation spectrum of 
noises can be reconstructed. This allows not only for 
checking for presence of a common component in the 
noises, but also for investigation of signal propagation 
through the environment, and of causal correlations be¬ 
tween noises affecting the two qubits. 
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Appendix A: Filtering term 
1. Dirac delta-like approximation to filtering term 

The PDD and CP sequences are defined as 


/p^dW = E(-i)'^n ^ = E(-i)'=n 


k=0 


fc=0 


n + 1 
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t — 
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kT 


2 {n -1-1) n + 1J 


(Al) 
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n+l f kT 
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2(n+l) 

T 


t — 


T 
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T { 2 nTl)T 

4 (n-|-l) 2 (n-|-l) 


(A2) 


Here, n(a;) is a box function that is 1 for — 1 < x < 1 and 0 otherwise, T is the total time and n is the (odd) number 
of pulses. The Fourier transforms the filtering functions can be easily computed 
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/cp^’^(‘^) = e /poo(w) - e*2(n4i) _:_sinc 
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Now let’s consider the first maximum of the modulus square of the filtering term |/pbD(^)l^ 


fa,T I _ {n + l) 7 r 
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(AS) 


The peak is very sharp for large n. This is because at u = (n + l)7r/r the sin function in the numerator oscillates 
n + l times faster then any other component of /. Consequently, this sine averages to ^ ^ and there is nothing left 
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to compensate the very small factor l/(n + 1)^ multiplying everything. This behavior repeats itself for each sector 

a; G ( (m — 1) (m + 1) ) with m = ±1, ±3, ±5,_ Hence, it is reasonable to approximate the ratio 

of sin and cos with delta function and simply set all other factors as constants at ujm = rn{n + 1 )tt/T. With this 
approximation the damping rate of the coherence is given by: 


Xaa = ^ y <S'aa(w)|/pijD(a;)P^ = J S'aa(w) 
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m=±(l,3,...) 

where we used the relation 5'(—tc) = S*{u!). By inspecting Eq. 


_ 4r J e 
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we can now identify A(z) function: 
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The imaginary part of cross spectrum is revealed when the filtering term is a product of both sequences: 

Xi2 = ReJ fpuui^) ‘^i2(w)^=y |/piy,(a;)pS'i2(a;)e*5I^E + o((,^ + l)-i) 
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2. Discussion of the approximation 

We estimate that the corrections to “delta-like” approximation scales like (n-|- 1)“2 m (^ = 1, 2,...) and are insignif¬ 
icant for large enough n. However, the same cannot be said about corrections due to additional peaks of sine function 
(the sum over m in Eq. (IA8D 1. Hence, if we want to precisely estimate the power spectrum S{uj) at some point lu it 
is not enough to simply read out the coherence at given time To and apply the proper number n of pulses. Since the 
damping rate is a combination of S-s at many different, successionally higher frequencies, the measurement has to 
be complemented by additional data points at shorter and shorter times [l6| . Here we adopt the following notation: 

Sk = S {2k -\- 1) and Xk = X ( 2 M -1 ’ coherence damping rate at time To/{2k + 1) with n pulses; then 


the measured decay rates and the values of the spectral density satisfy the following linear relation 

Xk = f^uif/^^sr 

1=0 

Here {S^) is the real (imaginary) part of the spectral density and the ^max x Z^ax matrix is defined as 
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The relation (|A 9 |) can be inverted, and as a result we obtain the approximation for real or imaginary part of spectral 
density 
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(A 13 ) 


For example, in our numerical experiment we set /max = 13 , that is, we assumed Siio > 27 x (n + l)7r/ro) ^ 0 . In 
such a case the explicit form of matrix jg given by 
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Appendix B: Correlation functions 


1. Coupled telegraph noises 


We define the system of coupled dichotomic fluctuators by listing all quantities W{a2,h2\ai^hi) - the probability 
per unit time of the transition 


{A in state oi AND B in state bi) —> {A in state 02 AND B in state 62), (Bl) 

W( 1 ,&| 0 ,&) = 7 A (B 2 ) 

W (0 , & 11 , &) = 7 A (B 3 ) 

W(a, l|a, 0 ) = 7 ft (^ 4 ) 

W(a, 0 |a,l) = 7 f^ (B 5 ) 


The probabilities for other process either vanish or are simply irrelevant. 
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2. Steady-state probability distribution 

The steady-state probability distribution is given by 

pTi = — [ 7 . 1 ( 7 ^ + 7n) + 7n(7^4. + "lu)] > 

[ 7 ^( 7 ?^ + 71®) + 7u(7®4, + 7®)] , 

Pfo = [ 7 ^( 7 ^ + 7i®) + 7^(7®t + 7®)] . 

Pm = [ 7 ^( 7 ®^ + 7i®) + 7®t(7®4, + 7®)] , 

where 

u+ = 27a (7(u + 7® + 7(^ + li-\) + 2(7(^ + 7^) hoi + 7u) ■ 


3. Cross-correlations 


The cross-correlation functions are given by fairly complicated formulas 
CABir > 0 ) = (a(r) 6 ( 0 )) = 
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4. Cross-power spectrum 


For brevity sake, let’s rewrite the correlation functions as follows 
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Now we wish to calculate the Fourier transform of those functions 
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where we used the relation CapiT > 0 ) = C^ai—r). The explicit form of the cross-power spectrum is given by 
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The above formulas simplify significantly when Ays = 0 : 
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Appendix C: Details of the numerical experiment 

The parameters of the coupled telegraph noises were set to: 7^ = 7^ = 0 , 7^ = 7^ = 4 , 7^ = 1 and ui = U2 = 0.1 . 
In such a case we have Ays = 0 and e = — 16 . The number of pulses n = 37 have been kept fixed and the total time 
T was manipulated in order to sweep the wide range of frequencies. 



